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ABSTRACT 

We report on the prompt X- and 7-ray observations of GRB990712 with the 
BeppoSAX Gamma-Ray Burst Monitor and Wide Field Camera No. 2. Due 
to Sun constraints, we could not perform a follow-up observation with the Bep- 
poSAX Narrow Field Instruments. The hght curve of the prompt emission shows 
two pulses and a total duration of about 40 s in X-rays. In gamma-rays the 
event is even shorter. The 2-700 keV spectral emission with time shows a dis- 
continuity in the peak energy Ep of the EF(E) spectrum: Ep is above our energy 
passband during the first pulse and goes down to ~ 10 keV during the second 
pulse. Another peculiarity is noted in this event for the first time: the evidence 
of a 2 s duration emission feature during the tail of the first pulse. The feature 
is consistent with either a Gaussian profile with centroid energy of 4.5 keV or 
a blackbody spectrum with kT^h ~ 1.3 keV. We discuss the possible origin of 
the feature. The most attractive possibility is that we are observing the thermal 
emission of a baryon-loaded expanding fireball, when it becomes optically thin. 

Subject headings: gamma rays: bursts — gamma rays: observations — X-rays: 
general — shock waves 

1. Introduction 

Observations of cosmic Gamma- Ray Bursts (GRBs) with the BeppoSAX satellite are 
providing a key contribution to theories about their nature. Among the still unsettled 
questions, it is still not clear what are the mechanisms that produce the observed X-ray 
spectra and their evolution with time (see, e.g., Frontera et al. 2000) and what are the 
environments in which the GRBs occur. In the context of the internal shock model. 
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synchrotron radiation is generally expected to play an important role in the production 
of the observed GRB spectra (e.g., Tavani 1996), but Inverse Compton can also give 
a significant contribution to them (Ghisellini et al. 2000). Also blackbody emission 
from the photosphere of the fireball (Meszaros & Rees 2000) is expected to contribute 
to the GRB spectra, and inhomogeneities in the GRB outfiow, made of dense highly 
ionized metal-rich material, could give rise to broad spectral features, mainly K-edges 
(Meszaros and Rees 1998). Effects of photoelectric absorption and Compton scattering 
from the circumburst material can modify the intrinsic energy spectrum of the GRBs, 
with the introduction of absorption cut-offs and features, such as K-edges and emission 
lines (Meszaros and Rees 1998; Bottcher et al. 1999), the presence of which has already 
been reported for some GRBs (Yoshida et al. 1999; Piro et al. 1999; AntoneUi et al. 2000; 
Piro et al. 2000; Amati et al. 2000). The separation of the intrinsic and external components 
is of key importance for establishing both the GRB emission mechanisms and the properties 
of the GRB environment. 

The Gamma Ray Burst Monitor (GRBM, 40-700 keV, Frontera et al. 1997), and 
the two Wide Field Cameras (WFC, 2-28 keV, Jager et al. 1997) on board BeppoSAX 
offer the opportunity to study the GRB energy spectra in the 2-700 keV energy band 
where the above components can be investigated (e.g., Frontera et al. 2000). GRB990712 
was detected by the WFC No. 2 and the GRBM, showing in the 2-26 keV band the 
highest peak fiux ever observed from a GRB with BeppoSAX. Its position was promptly 
distributed to the astronomical community (Frontera et al. 1999). Follow-up observations 
with the BeppoSAX Narrow Field Instruments were not possible due to Sun constraints. 
Observations were performed in the optical and radio bands. An optical transient (OT) with 
magnitude R — 19.4 ± 0.1 was discovered about 3 hrs after the event (Bakos et al. 1999) 
and its redshift is now well determined {z — 0.4331 ± 0.0004) (Vreeswijk et al. 2001). 
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2. Observations 

GRB990712 was detected on July 12, 1999 starting on 16:43:02 UT (Prontera et al. 1999). 
Its position was determined with an error radius of 2' (99% confidence level) and was 
centered at ctsooo = 22^31'^50^ ^2000 = -73°24'24" (Heise et al. 1999). Features and 
data available from GRBM and WFCs have already reported in several papers (e.g., 
Prontera et al. 2000 ). The effective area exposed to the GRB was ^ 420 cm^ in the 
40-700 keV band and 37 cm^ in the 2-26 keV energy band. The background in the WFC 
and GRBM energy bands was fairly stable during the event. The WFC spectra were 
extracted through the Iterative Removal Of Sources procedure (IROS, e.g. Jager et al. 1997 
) which implicitly subtracts the contribution of the background and of other point sources 
in the field of view. The count rate spectra were analyzed using the XSPEC v. 10 software 
package (Arnaud 1996). The quoted errors for the spectral parameters correspond to 90% 
confidence. Parameter values shown in brackets in Table 1 have been fixed while fitting. 



3. Results 

Figure 1 shows the time profile of GRB990712 in four energy bands after the 
background subtraction. In all bands the GRB shows a double-pulse structure, with an 
opposite behavior with energy of the first pulse with respect to the second: the peak fiux 
of the first pulse increases with energy, while that of the second pulse decreases. The total 
duration of the event is about 20 s above 100 keV, but much longer (about 40 s) over the 
full 2-26 keV range. The spectral evolution of the event was studied by subdividing the 
GRB time profile into 8 adjacent time intervals and performing an analysis of the spectrum 
of each interval (see Fig. 1). We fit the spectra with either a power law (PL, N{E) oc E"), 
a broken power-law (bkn-pl) or a smoothly broken power law (bl. Band et al. 1993), 
with photoelectric absorption (wABS model in XSPEC) . We assumed the Galactic hydrogen 
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column density A^h (= 4.52 x 10^° cm"^, Dickey & Lockman 1990 ) along the line of sight 
to the GRB. For the first three time shces the fit was also performed leaving A^h free to 
vary, but unfortunately A^h was not constrained by the data. The fit results with a PL 
and a BL are given in Table 1. The results obtained with the bkn-pl model were similar 
to those obtained with the BL model and are not reported in Table 1. The spectra of the 
time intervals A and B, that correspond to the rise of the first pulse of the GRB, and the 
spectrum of the interval D, that corresponds to the late tail of this pulse, are well fitted 
with a PL model. The BL model provides a good description of the spectra in the time 
intervals E and F, that correspond to the core of the second pulse, while the PL model 
again well describes the tail of the second pulse (intervals G and H). However neither of 
these models provides a good description of the C spectrum {y^ jv — 16/6 for a PL and 
16/5 for a bl). An excess with a significance level of about 3(7 is evident around 4 keV. We 
investigated a possible instrumental origin of this anomaly with negative results. The WFC 
high-voltage, which is monitored every second, does not show any glitch in the interval C. 
Given the accuracy of the readout (which dominates over statistical noise) this excludes 
gain changes higher than 0.01%. As in other GRB detections, the WFC ratemeters all 
show the gamma-ray burst, even the ones that measure the illegal events, so we do not 
find any anomaly in this event. There are no dips or spikes of any sort, with a typical 3(7 
upper limit of 10% per second of measurement (for a count rate of 700 counts s~^). There 
is one thermometer that shows a small change about 10 s before the burst, but all other 
measurements do not show anything out of the ordinary. The count excess in the GRB 
spectrum of the interval C is also clearly visible (see Figure 2) in the ratio between the 
C count spectrum and the Crab spectrum measured when this source was observed at an 
angular offset similar to that of GRB990712. For comparison. Figure 2 also shows the ratio 
with Crab of the spectra measured in the intervals B and D, that precede and follow C, 
respectively. As can be seen, in the interval B the greater hardness of the GRB spectrum 
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( photon index of ~1.4, see Table 1) than that of Crab is apparent, while in the interval D 
the flatness of the Crab ratio is consistent with the similar slope of the GRB spectrum (see 
Table 1) with that of Crab. A slight hint of the 4 keV excess in the spectrum also appears 
in the interval B, but it is not statistically significant. We point out that the Crab ratio 
technique is adopted, to discover cyclotron lines in the spectra of X-ray pulsars (e.g., Dal 
Fiume et al. 2000). 

The addition of a Gaussian function or a blackbody spectrum (bb) to a PL model 
provides a good fit ix^ /v = 1.6/5 and x^/z^ = 6.2/6, respectively) of the C spectrum. The 
best fit parameters of both models are reported in Table 2. For a better determination of 
the Gaussian and BB parameters, the photon index a of the PL model was kept fixed in 
the fit to the best fit value, that is given by 1.34 ± 0.17 or 1.24 ± 0.20, depending on the 
model assumed for the feature, a Gaussian or a BE model, respectively. The count rate 
spectrum of the interval C along with the best fit curve of the PL plus a blackbody model 
is shown in the top panel of the Figure 3, while the ratio between the count spectrum and 
the best fit power-law alone is shown in the bottom panel. The excess counts to the PL 
model are apparent. The evolution of the logarithmic power per photon decade (the EF{E) 
spectrum) with the time from the GRB onset is shown in Figure 3. The emission feature 
during the C interval is also apparent in this plot. The peak energy of the EF{E) spectrum 
(see Table 1) is above our energy passband for the entire duration of the first pulse before 
suddenly becoming much lower (~ 10 keV) from the beginning of the second pulse. 

From the spectral fits we derived GRB fiuence and peak fiux. The 7-ray (40-700 keV) 
fiuence of the burst is = (6.5±0.3) x 10~^ erg cm~^, while the corresponding value found 
in the 2-10 keV band is Sx = (2.60±0.06) x 10"^ erg cm-^, with a ratio Sx/S^ = 0.40±0.03, 
which is one of the highest values found with BeppoSAX (Frontera et al. 2000). The 
2-700 keV fiuence is given by = (1.10 ± 0.03) x 10"^ erg cm'^. The 7-ray peak fiux 
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is P-y = 4.1 ± 0.3 photons/cm^ s corresponding to (1.3 ± 0.1) x 10 ^ erg cm ^ s ^, while 
the corresponding 2-10 keV peak flux is Px = 41 ± 4 photons/cm^ s, corresponding to 
(3.3 ± 0.3) X 10"^ erg cm'^ s'K 



4. Discussion 

From the redshift value of the optical afterglow of GRB990712 {z — 0.4331, Vreeswijk 
et al. 2001) we can derive the X- plus 7-ray energy released in the main event. Assuming 
isotropic emission and a standard Friedman cosmology {Hq = 70 km s~^ Mpc~^ and go = 0.5), 
we get a 2-700 keV released energy of Ej-ei = (5.9 ± 0.2) x 10^^ erg. A sizeable fraction 
of this energy (~20%) is released between 2 and 10 keV. If we exclude the controversial 
case of GRB980425/SN1998bw (Galama et al. 1998; Kulkarni et al. 1998; Pian et al. 2000), 
GRB990712, in addition to showing the lowest redshift, is one of the least energetic events. 

GRB990712 is marked by a peculiarity, which is noted here for the first time: the 
evidence of a broad emission feature at 4.4 keV, visible for 2 s, superposed on a power 
law continuum model. It can be described by either a Gaussian profile with a full width 
at half maximum of ~3 keV or a blackbody emission with kT^f, — 1.3 keV. If we assume 
the Gaussian description, its centroid energy corrected for redshift {Eq = 6.4 ±1.1 keV) is 
consistent with the energy of both an iron K fiuorescence line and a iron recombination line. 
The interpretation of the emission feature as an iron recombination line is tempting, yet 
it makes some stringent demands on models: in fact, it requires much mass to be present 
within a few light seconds of the burst site, and for this mass to be moving at Newtonian 
speeds. Vietri et al (2001) derive the expected rate of photons for a narrow line: 

4x10^2 Mp^ 6x10^^ cm 

where T7 and R are the electron temperature (units of 10^ K) and the external radius of the 
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line emitting medium, respectively, and Mpe is the iron mass present in this medium. For 
a broad hne, such as that in GRB990712, the above value is an understimate by a factor 
of a few at most. Assuming Ty = 1, i? = 6 x 10^^ cm and Mpe = IMq, comparison with 
Table 2 shows that Eq. 1 underestimates the observed hne luminosity by four to five orders 
of magnitude; inserting R — 10 light seconds in the above equation yields the correct line 
luminosity, but for a total mass of iron of at least Mp^. = O.IMq. Assuming a realistic iron 
relative abundance (at least 1% of the total mass, if what we are seeing is a type I SN, but 
more for any other hypothesis) shows that we must explain the presence of at least IOMq 
of matter, at radii of a few light seconds, with none of this obstructing the line of sight. 

A more attractive possibility is that the observed feature is indeed thermal. Though we 
cannot definitely establish the thermal nature of the emission during the C interval, we wish 
to remark that the fireball model can account naturally for the presence of these features. 
In fact, as remarked by both Paczynski (1986) and Goodman (1986) hyper-relativistic 
expansion naturally leads to the liberation of a fully thermal spectrum, at the time when 
the fireball expansion becomes optically thin. The initial absence, and later disappearance, 
of the peak in question does not create difficulties within the fireball model: it can easily be 
ascribed to inhomogeneities in the time-structure of the relativistic wind, inhomogeneities 
which are in any case required in order to account for the burst sub-second variability. 
The fact that the peak, furthermore, appears during the tail of the first pulse of course 
makes the detection of the weaker thermal component easier (see the very revealing Fig. 
2 of Meszaros and Rees 2000). If this spectral feature is indeed thermal in origin, within 
the fireball model there is quantitative, and independent check on this identification. In 
fact, the expected photospheric radius within the fireball model is (Meszaros & Rees 2000) 
Vph — 1.2 X 10^^ cm L52Yrj2^, where L52 is the wind luminosity in units of 10^^ erg s"^, 
y ~ 1 is the number of electrons per baryon, and ri2 is the flow Lorentz factor rj in 
units of 100. At this radius, the observed photospheric luminosity and temperature are 
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Lph = L52{v/v*)^^^ and Qph = Qo{r]/ri^f/^, respectively, where 77* ^ 10^(L52//i ^yFo)^/^ and 
/^i Tq Y''^ keV, Fq (> 1) being the initial bulk Lorentz factor of the wind, 
and III the mass, in units of 10 Mq, of the rotating black hole, from 6 times the gravitational 
radius of which the fireball is assumed to start its expansion. By fitting simultaneously 
the observed temperature corrected for the redshift (1.86 keV) and the luminosity of the 
photosphere (~ 2 x lO^o erg s-^), we find r] ^ 100Y^/^{fi^^Tl^y/''^ and L52 ^ 2{fiiTl^)'^/^ . 
We thus see that the two independently determined observational parameters, blackbody 
temperature and luminosity, are well-fitted by values of the theoretical parameters, 77 and 
L52, well within the expected range. Assuming unit values for Y, fii and Fq, L52 results 
to be about 100 times higher than the estimated 2-700 keV luminosity (~ 2 x 10^° erg 
s~^) assuming isotropy. That would imply an efficiency of only 1% in the production of 
electromagnetic radiation. 

In addition to the transient feature, the event shows a spectral evolution that is not 
typical of other GRBs observed with BeppoSAX (Prontera et al. 2000): the peak energy Ep 
of the EF{E) spectrum (see Table 1 and Figure 3) is constantly above our energy passband 
for the entire duration of the first pulse, while it takes a low value (~ 10 keV) with the 
onset of the second pulse. This discontinuity can be the result of two successive electron 
acceleration episodes, giving rise to the first and the second pulses. The different peak fiux 
behavior of the two pulses with energy, discussed in section 3 confirms this scenario. The 
emission feature is found only during the first acceleration event. 

Thanks to John Stephen for his careful reading of the manuscript. Also many thanks 
to the anonymous referee who greatly stimulated us to improve the paper. 
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Fig. 1. — Light curves of GRB990712 in four energy bands, after background subtraction. 
The zero abscissa corresponds to 1999 July 12, 16:43:01.6 UT. The time intervals over which 
the spectral analysis was performed are indicated with vertical dotted lines. 

Fig. 2. — Ratio with the Crab count spectrum of the GRB990712 spectra in the time intervals 
B, C and D, respectively, measured with the WFC No. 2. The Crab spectrum used was 
measured when this source was observed at an angular offset similar to that of GRB990712. 

Fig. 3.— WFC + GRBM spectrum of GRB990712 in the time interval C, along with the 

best fit curve obtained assuming a power-law plus a blackbody model. Bottom panel: ratio 
between the data and the best fit power-law model alone. 

Fig. 4. — EF[E) spectrum of GRB990712 in the time intervals in which we divided the burst 
time profile (see also Table 1). It is apparent the GRB pecuhar spectrum in the interval C. 
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Table 1. Spectral evolution of GRB990712 prompt emission 



Slice 


Duration (s) 


Model^^: 




a 




Ep(keV) 




A 


4 


PL 


3.2±2.5 


-1.40 ±0.09 




> 700 


7.8/6 






PL 


[0.0452] 


-1.34 ±0.07 






10./7 


B 


2 


PL 


2.7±2.1 


-1.44 ±0.08 




> 700 


4.4/5 






PL 


[0.0452] 


-1.38 ±0.06 






7/6 


C 


2 


PL 


0.3±2.0 


-1.66 ±0.11 




> 700 


16./6 






PL 


[0.0452] 


-1.64 ±0.07 






16/7 


T— \ 

D 


4 


PL 


[0.0452] 


—1.80 ± 0.07 




> 700 


4.6/6 


E 


2 


PL 


[0.0452] 


-2.04 ±0.05 






37/7 




2 


BL 


[0.0452] 


-0.4 ±0.7 


-2.4 ±0.3 


11±8 


3.5/5 


F 


3 


PL 


[0.0452] 


-2.01 ± 0.04 






89/7 




3 


BL 


[0.0452] 


-0.7 ±0.4 


-2.3 ±0.2 


7±3 


4.3/5 


G 


7 


PL 


[0.0452] 


-2.15 ±0.05 






11/7 




7 


BL 


[0.0452] 


-1.6 ±0.6 


-2.3 ±0.3 


20±15 


2.5/5 


H 


7 


PL 


[0.0452] 


-2.3 ±0.2 






3.0/7 



("^The BL (Band Law) refers to the smoothed broken power-law proposed by Band 
et al. (1993): a and /3 are the power-law photon indices below and above the break 
energy Eq, respectively, and Ep=Eo(2±ci;) is the peak energy of the EF(E) spectrum. 

(''^A^H values are given in units of 10^^ cm~^. 
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Table 2. Best fit parameters of tlie time slice C spectrum 



Parameter 


PL+Gaussian 


PL±BB 


a 


[1.34] 


[1.24] 


FpL(@lkeV) [cm-2 g-i] 


1.95 ±0.24 


1.13 ±0.14 


£^iine [keV] 


4.4 ±0.8 




CTiine [keV] 


1.4 ±0.7 




/line [10-^ erg cm-2 s"^] 


2.7 ± 1.1 




Lii,e [10^^ pilot S-1] 


2.5 ±0.9 




A;rbb [keV] 




1.3 ±0.3 


[10^9 erg s-i] 




2.5 ±0.6 




1.6/5 


6.2/6 
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Fig. 4 



